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Abstract 

The spectral analysis of the electromagnetic field on the background of a infinitely thin flat 
plasma layer is carried out. This model is loosely imitating a single base plane from graphite and it 
is of interest for theoretical studies of fullerenes. By making use of the Hertz potentials the solutions 
to Maxwell equations with the appropriate matching conditions at the plasma layer are derived and 
on this basis the spectrum of electromagnetic oscillations is determined. The model is naturally 
split into the TE-sector and TM-sector. Both the sectors have positive continuous spectra, but the 
TM-modes have in addition a bound state, namely, the surface plasmon. This analysis relies on the 
consideration of the scattering problem in the TE- and TM-sectors. The spectral zeta function and 
integrated heat kernel are constructed for different branches of the spectrum in an explicit form. 
As a preliminary, the rigorous procedure of integration over the continuous spectra is formulated by 
introducing the spectral density in terms of the scattering phase shifts. The asymptotic expansion 
of the integrated heat kernel at small values of the evolution parameter is derived. By making 
use of the technique of integral equations, developed earlier by the same authors, the local heat 
kernel (Green's function or fundamental solution) is constructed also. As a by-product, a new 
method is demonstrated for deriving the fundamental solution to the heat conduction equation (or 
to the Schrodinger equation) on an infinite line with the 5-like source. In particular, for the heat 
conduction equation on an infinite line with the (^-source a nontrivial counterpart is found, namely, 
a spectral problem with point interaction, that possesses the same integrated heat kernel while the 
local heat kernels (fundamental solutions) in these spectral problems are different. 
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I. INTRODUCTION 



Now the employment of the fuUerenes in practical activity becomes of rising impor- 
tance (in superconducting technique, in mechanical engineering for reducing the friction, in 
medicine for drug production and so on). FuUerenes 0] are the third crystal modifica- 
tion of carbon. First they were produced in 1980-ties (Nobel price in Chemistry in 1996, 
R. F. Curl, H. W. Kroto, and R. E. Smalley). The elementary blocks of the fuUerenes are 
the giant carbon molecules (Ceo, C70, C84, and greater) having the form of empty balls, el- 
lipsoids, tubes and so on. There appears a special branch of technology, nanotechnology , 
which is based on using the fuUerenes and especially, nanotubes. In theoretical studies of 
the fuUerenes it is important to estimate the ground state energy (the Casimir energy) of 
such molecules. The carbon shells, forming the structure of fuUerenes, can be considered as 
infinitely thin layers, and their collective electrons can be treated as two-dimensional plasma 
layer by making use of the hydrodynamic picture for its description. As a result, the prob- 
lem reduces to the solution of the Maxwell equations with charges and currents distributed 
along the surfaces modeling the carbon sheets. The method of boundary (matching) condi- 
tions applies here, i.e., the Maxwell equations are considered only outside the plasma layers 
and when approaching these sheets the electric and magnetic fields should meet appropriate 
matching conditions. From the mathematical stand point the pertinent spectral problem 
proves to be well posed. 

For constructing the quantum version of this model and, first of all, for a rigorous treat- 
ment of the divergencies encountered here the spectral analysis of the underlying field theory 
should be accomplished, i.e., the spectrum of the dynamical system in question should be 
determined, the mathematically consistent integration procedure over this spectrum should 
be defined and on this basis the spectral zeta function and heat kernels (integrated and local 
ones) should be calculated as weU as their asymptotic expansions under small values of the 
evolution parameter should be found. 

The present work is exactly devoted to the spectral analysis of an infinitely thin plasma 
sheet model with the simplest geometry, namely, the plasma layer has the form of a plane. 
The spectrum of the model proves to be very reach, namely, it contains both continuous 
branches and bound states (surface plasmon), and we managed to construct all the pertinent 
spectral functions possessing interesting properties. 

The outline of the paper is as follows. In Sec. [H] the physical formulation of the model 
under study is briefly given. In Sec. Illll the solutions to the Maxwell equations obeying the 
pertinent matching conditions at the plasma layer are constructed. The use of the Hertz 
vectors for this aim proves to be very effective. Proceeding from this the spectrum of elec- 
tromagnetic oscillations in the model is found. As usual in electrodynamical problems the 
model is naturally separated into the TE-sector and TM-sector. Both the sectors have pos- 
itive continuous spectra, but the TM-modes have in addition a bound state, namely, the 
surface plasmon. This analysis relies on the consideration of the scattering problem in the 
TE- and TM-sectors. In Sec. IIVI the spectral zeta function and integrated heat kernel are 
constructed for different branches of the spectrum in an explicit form. As a preliminary, the 
rigorous procedure of integration over the continuous spectra is formulated by introducing 
the spectral density (the function of the spectral shift) in terms of the scattering phase shifts. 
Sec. 13 is devoted to the construction of the local heat kernel (Green's function or fundamen- 
tal solution) in the model at hand. Here the technique of integral equations governing the 
heat kernel to be found is used. As a by-product, a new method is demonstrated here for 
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deriving the fundamental solution to the heat conduction equation (or to the Schrodinger 
equation) on an infinite line with the S-like source. In particular, for the heat conduction 
equation on an infinite line with the 5-source a nontrivial counterpart is found, namely, 
a spectral problem with point interaction, that possesses the same integrated heat kernel. 
However the local heat kernels in these spectral problems are different. In the Conclusion 
(Sec. IVI|) the obtained results are summarized briefly. 



II. FORMULATION OF THE MODEL 

In recent papers 0, 0] Barton has proposed and investigated the model of infinitesimally 
thin two-dimensional plasma layer that is loosely inspired by considering a single base plane 
from graphite or the giant carbon molecule Ceo- Effectively the model is described by the 
Maxwell equations with charges and currents distributed along the surface S : 

V-B = 0, V X E-icuB/c = 0, (2.1) 

V ■ E = 47r5(x - xe) a, V x B + icuE/c = 47r5(x - xs) J/c . (2.2) 

It is assumed that the time variation of all the dynamical variables is described by a common 
factor e"*"^*. 

The properties of the plasma layer (for example, plasma oscillations and screening) differ 
considerably from those in bulk media In hydrodynamic approach the plasma is con- 
sidered as the electron fluid embedded in a rigid uniform positive background. Electrical 
neutrality requires that the equilibrium electron charge density —erio precisely cancels that 
of the background. In a dynamical situation, the electron density is altered to r2(x, t), and 
the equation of continuity requires 

^ + V-(nO = 0, (2.3) 

where ^(x, t) is the electron velocity at the point x. Usually one considers the /mear response 
of an initially stationary system to an applied perturbation. In this case, the induced charge 
density 

or(x,t) = -e(n(x,t) -no) , (2.4) 

the velocity and the fields E and B are of first order. In this approximation the continuity 
equation (j2.3p gives for the plasma layer 



cr 



e?2oV|| ■ ^ = or (j = enoV||-^. (2.5) 



The superscripts || and ± indicate here and below the vector components respectively parallel 
and normal to the surface S. 

For the induced charge current J we have 

3 = —erio^ = i eu no$, . (2.6) 

The Newton second law applied to an individual electron gives 

m ^'(x, t) = -e E|| (x, t), x G S, or ^ = En . (2.7) 
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Finally, the induced-on-the-surface charges and electric currents are determined by the par- 
allel components of electric filed 



,= £!!!|V„.E„ J^i'^E,. (2.8) 

mcj^ muj 

As known 0], electromagnetic field, generated by charges and currents with singular 
densities is described mathematically in the following way: outside the singularities the 
Maxwell equations (j2.ip . (|2.2p without sources should be satisfied 

V-B = 0, V X E-zc.;B/c = 0, (2.9) 
V-E = 0, V X B + ituE/c = 0, x^S, (2.10) 

and when approaching at singularities the limiting values of the fields should meet the 
matching conditions determined by the singular sources. In the regions of singularities the 
Maxwell equations are not considered. The singular charges and currents in right-hand sides 
of the Maxwell equations ()2.2j) lead to the following matching conditions 

[E||] = 0, [EJ = 2g(c/a;)2V|| • E|| , (2.11) 
[B±] = 0, [B||] = -2iq{c/uj) n x Ey . (2.12) 

Here g is a characteristic wave number q = 27me^/mc^, the square brackets [F] denote the 
discontinuity of the field F when crossing the surface S, and n is a unit normal to this 
surface usually used in formulation of matching conditions 0. 

The physical origin of the model leads to a Debye-type cutoff K on the surface-parallel 
wavenumbers of waves that the fluid can support. However if one digresses from this cutoff 
then the Maxwell equations (|2.9|) and (j2.1Up with matching conditions (|2.1H) and (|2.12j) can 
be regarded as a local quantum field model that is interesting by himself. First of all the 
question arises here concerning the analysis of divergences. As known [1, H,|l3, 11. Il^ ll3l|. 



to this end one has to calculate the heat kernel coefficients for the differential operator 
garnering the dynamics in this model (see Eqs. ()2.9|) - ()2.12j) ). 



III. SOLUTION TO THE MAXWELL EQUATIONS FOR A FLAT PLASMA 
SHEET 

When constructing the normal modes in this problem it is convenient to use the Hertz 
potentials 0, Q, • In view of high symmetry of the fiat plasma layer the electric (E) and 
magnetic (B) fields are expressed in terms of electric (11') and magnetic (H") Hertz vectors 
possessing only one nonzero component 

W = e,e'^^<^{z), n" = e,e^^^'^{z). (3.1) 

Here the plasma layer is taken as the [x, y) coordinate plane, the axes z is normal to this 
plane, k is a two-component wave vector parallel to the plasma sheet and s = (x, y); e^., ej,, 
and are unit base vectors in this coordinate system. The common time- dependent factor 
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For given Hertz vectors 11' and 11" the electric and magnetic fields are constructed by 
the formulas 

E = VxVxn', B = -i'^Vxn' (TM - modes,5^ = 0); (3.2 

E = z - V X n", B = V X V X n" (TE - modes, = 0) . (3.3 
c 

Substituting here 11' and 11" by Eq. ()3.1|) one obtains for the TE-modes 



E = i-kye, + Key) - e^'^^^(z), (3.4 



c 

B = i{k^e^ + kyey)e^^^^'{z) + e,k'^(^^^^{z) (3.5 



and for the TM-modes 



E = i{k^e^ + kyey)e^^^^'{z) + e,k'^(^^^^{z) (3.6 
B = (A;,e,-A;,e,)-e^^^$(z), (3.7 



The fields ()3.4|1 - ()3.7p will obey the Maxwell equations ()2.9j) and ()2.1()|l outside the plasma 
layer [z = 0) when the functions e'^^^^{z) and e^^^'^i^z) are the eigenfunctions of the operator 
(—A) with the eigenvalues co'^/c^ (A is the three-dimensional Laplace operator): 



We deduce from here 



^2 

-Ae^'^z) = —e'^'^z), (3.8 
-Ae^^^^(z) = — e^'^"^(z). (3.9 



^"{z) = ["^ - e]^z), (3.10 
^"{z) = ('4-k']n^), (3.11 



— oo <2<oo, -27^ 0, 0<A;^<oo. 

The substitution of the fields ()3.4|) - ()3.7|) into the matching conditions ()2.11|) and ()2.12|) 
gives the analogous conditions at the point z = for the functions ^{z) and \E'(2;): 

[$(0)] = -2g(^)V(0), [$'(0)] = 0, (3.12 
[vl/(0)] = 0, [vl/'(0)] = 2gvl/(0) , (3.13 

where the notation is introduced 

[F{z)]=F{z + 0)-F{z-0). (3.14) 

The matching conditions ()3.12|) are very specific, they involve the eigenvalues of the whole 
initial spectral problem ()3.8|) . It implies, specifically, that the dynamics of TM-modes along 
the z axes 'feels' the parallel dimensions {x, y). Thus for the TM-modes the initial spectral 
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problem does not factorizes completely into two independent boundary value problems along 
the z axis and in parallel directions. 

In order to complete the formulation of the boundary spectral problems ()3.10|) - ()3.13|) 
we have to specify the behaviour of the functions $(^) and ^[^z^ when \z\ oo. Proceeding 
from the physical content of the problem in question, we shall consider such functions $(-2) 
and ^1/(2;) which either oscillate (p^ = cj^/c^ — A;^ > 0) or go down (k^ = k"^ — uj^ j(? > 0) 
when — i> 00. In the first case we are dealing with the scattering states and in the second 
one the solutions to the Maxwell equations describe the surface plasmon.^ 

The spectral problem for the TE-modes ()3.1H) and ()3.13|) with the regarding boundary 
conditions is a selfadjoint spectral problem. Substituting the TE-fields fl3.4|) and ()3.5p into 
the initial Maxwell equations with 5-sources ()2.2p we arrive at the equation for "^{z) with 
the 5-potential (instead of Eq. 1)3.111) ) 



-^"{z) + 2q6{z)^{z) =p'^{z), — = k'+p\ (3.15) 



c 



Integration of this equation in vicinity of the point z = leads to the matching conditions 
(|3.13p (see, for example, 2l||). The potential in Eq. ()3.15|) is equal to 2q6{z). As it was 



mentioned above, the spectral problem for the TM-modes ()3.10p . ()3.12|) possesses a pecu- 
liarity, namely, the matching conditions ()3.12|) involve the eigenvalue u/c. Obviously, the 
selfadjointness condition is not satisfied here. 

The scattering problem for the functions ^{z) and "^{z) is formulated on a line: 
—oo<z<oo. As known [i^, 2^, the one-dimensional scattering problem has features 



in comparison with description of the scattering processes in three-dimensional space. 
The scattering states for the TE-modes are described by the functions 

^(z) = Cie'P' + C2e-'P', z<0, (3.16) 
^(z) = C3e'P', z>0, p>0, = — -k\ (3.17) 

The initial wave is coming from the left. Matching conditions ()3.13p give the following 
relation between the constants Q {i = 1, 2, 3): 

Ci + C2 = Cs, (3.18) 
^p{Cs-C^ + C2) = 2qCs. (3.19) 

From here we obtain for the reflection and transmission coefficients 

^TE _ ^ ^ ^ ^gi^^g.,^ ^TE ^ ^ cosr/e*\ (3.20) 

Ci p + iq p + iq 

where 77 (p) is the phase shift for the scattering of the TE-modes. When initial waves is 
coming from the right we obtain the same reflection and transition coefficients. There is no 



^ In electrodynamics of continuous media the surface plasmon is a special solution to the Maxwell equations 
which describes the electric and magnetic fields exponentially decreasing in the direction normal to the 
interface between two material media [i^IitIIi^ . The plasmon solution exists only in the TM-sector. It is 
these configurations of electromagnetic field that contribute to the Casimir force between two semi-infinite 
material media 0, . 
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other conditions on p besides the positivity p > (continuous spectrum). The scattering 
matrix, S{p), is determined by the phase shift rj^p) 

tanr]{p) = -^, p>0 (3.21) 

in a standard way 

S{p) = e^* (3.22) 

Between TE-modes there is no solutions which go down when \z\ oo. Indeed, such 
solutions, in accord with Eq. p.llj) would have the form 



00^ 



^i/{z) = C2e-^', z>0, K = +^k^- — >0. (3.23) 

Substitution of Eq. ()3.23|) into the matching conditions ()3.13p gives 

Ci = C2 and n = -q < . (3.24) 

Thus, there is no solutions of the form ()3.23j) with positive n. This is in accord with the 
known fact [IH that the potential defined by the Dirac (5-function \S{z) leads to bound 
state only for A < 0. In our case A = 2g > (see Eq. (j3.15|) ). 
Finally, the TE-modes have the spectrum 



w (k,p) 



C2 



k^+p% 0<p<oo. (3.25) 



Here the contribution k"^ is given by the free waves propagating in directions parallel to 
the plasma layer, and the contribution p"^ corresponds to the one- dimensional scattering in 
normal direction with the phase shift ()3.2H) . 

Proceeding in the same way one obtains for the scattering of the TM-modes: 

n^^ = ^-p-— = -^sin^e^^, T™ = ^^^-t^ = cos e- , (3.26) 
p'' + k"^ + tpq p'' + k'' + tpq 

where the phase shift fi{p, k) is defined by 

tanMp,A;) = -^^, p>0. (3.27) 

The scattering matrix is defined here by the formula analogous to Eq. (j3.22|) : S{p,k) = 

It turns out that between TM-modes there are solutions which not only oscillate when 
\z\ — *• oo, but also go down in this limit. In this case the function ^{z) is defined by the 
equation 



$"(^) -(fc^-^ )$(;.)= 0, (3.28) 

where 



C2 



^2 

P - — = k'^ >0. (3.29) 



The solution we are interested in should have the form 



$(^) = Cie"", z<0, 

$(2) = C2e-'^^ z>0, (3.30) 



2 

sp 

,2 



where k = +a/A;2 — lu'^/c^ > 0. The matching conditions (j3.1'2j) give rise to the equation 
for k: 

K^ + Kq-k'^ = 0, (3.31) 

the positive root of which is 

«:=y^7p-|. (3.32) 
For the respective frequency squared we derive by making use of Eq. (|3.29|) 

^ (v7+4P-g) >0. (3.33) 

Thus, the frequency of the surface plasmon is real, and the solution obtained oscillates in 
time instead of damping. 

Finally, the spectrum of the TM-modes in the problem under study have two branches: i) 
the photon branch and ii) the surface plasmon branch. The photon spectrum of TM-modes 
is defined by 

^i^^ = k'+p', keM^ 0<j9<oo. (3.34) 

Here the contribution k^, as in Eq. ()3.25p is given by the free waves propagating in directions 
parallel to the plasma layer and p"^ is the contribution of the one-dimensional scattering with 
the phase shift ()3.27|) in normal, to the plasma sheet, direction. The surface plasmon branch 
of the spectrum is described by 



^sl(k) 



C2 



I (v/g2 + 4F-g) , keR\ (3.35) 



In the spectral problem ()3.10p . ()3.12|) with negative q the parameter k defined by Eq. 
()3.32|) is positive, however the respective frequency squared is negative 



a;2 



^ = -^(g' + klv/?+4P). (3.36) 
Hence, in this case we have a resonance state instead of the bound state (surface plasmon). 



IV. SPECTRAL FUNCTIONS IN A FLAT PLASMA SHEET MODEL 

In the previous section we have determined the spectrum in the model under consideration 
(see Eqs. fl3.25|) . ()3.34|) . and ()3.35p ) therefore we can proceed to construction of the spectral 
functions in this model, namely, the spectral zeta function 

as)=TrL-^ = J2K' (4.1) 



8 



and (integrated) heat kernel 



K{t) = Tr (e-*^) = ^ e"^"* , (4.2) 

n 

where L is the differential operator in Eqs. ()3.10|) and ()3.11|) with matching conditions at 
z = ()3.12|) . ()3.13p and conditions specified above for — ^> oo; A„ are the eigenvalues given 
by Eqs. ()3.25|) . ()3.34p . and ()3.35p . In Eq. ()4.2|) t is an auxiliary variable, < t < oo. In our 
case L = —A and t has the dimension [length]^. 

As a preliminary, it is worth stipulating the procedure of spectral summation entered 
Eqs. (|4.ip and (|4.2p . Obviously the contribution to the spectral functions generated by free 
waves is given by a simple integration over the respective wave vector k: 

d^k rkdk 

(4.3) 



(27r)2 J, 271 

The contribution of the scattering states we define by making use of the phase shift 
method, namely, the integration over dp will be carried out with the spectral density p{p) 
(the function of spectral shift): 

where 6{p) is the phase shits for TE and TM modes given by Eqs. ()3.21|) and ()3.27p . respec- 
tively. 



A. TE modes 

Taking into account all this we get for the spectral zeta function in the TE-sector of the 
model 



(27r)2 '^^ 71 dp ""^^ 27T^ Jo p^ + q^Jo {k^ + p^Y ' 

The left-hand side of Eq. (14. 5p is defined in the semi-plane Re s > 1 and this representation 
can be analytically extended all over the complex plane s save for separate points. In order 
to do this it is sufficient to express the integrals in Eq. ()4.5|) in terms of the gamma-function 

.TE. ^ q'-'' T{3/2-s)T{s-1/2)T{s-1) g^-^^ 

Stt"' T[s) 87r(l — s) cos(7rs) 

The integrated heat kernel for TE-modes is given by 

K^^it) = K^^=^\t) ■ K^'^=^\t) , (4.7) 

where K^^~'^^ (t) is the free two-dimensional heat kernel responsible for free waves propagating 
in directions parallel to the plasma layer 
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and 



71 
1 

r 



dpe 



,2+ d 
dp 



arctan — 



P 



p2 _|_ 



dp 



1 — erf ( q\/t 



i e*''^ erfc {qVt 



(4.9) 



where erf [qVt] is the probabihty integral 24. l25l| . 

In the plasma sheet model under consideration the parameter q is strictly positive. How- 
ever the spectr al p roblem can be considered for the negative q also when it possesses 
a bound state [2l| . The contribution of this state should be taken into account when calcu- 
lating the respective heat kernel i^q<o(t). But this heat kernel can be obtained by analytical 
continuation of Eq. ()4.9|) 



Kq<oit) = -e"^ [l-erf(^-|g|v^ 



- e 
2 



1 + erf (\q\Vt 



(4.10) 



Comparison Eqs. ()4.9|) and ()4.10|) gives the following simple relation between the integrated 
heat kernel for 6 potentials with positive and negative coupling constants 



(4.11) 



where Kgyoit) is given by Eq. ()4.9|1 . 

The structure of divergences in quantum field theory is determined by the coefficients of 
the asymptotic expansion of the heat kernel when t — > +0 



K{t) = {Ant)-'^/'- ^"^'^^ 



n/2 



ES, 



(4.12) 



71=0,1,2,. 



where d is the dimension of the configuration space in the problem at hand, ES stands for 
the exponentially small corrections. The coefficients Bn/2 for TE-modes can be derived in 
two different ways: i) by making use of Eqs. ()4.7|) - ()4.9|) which determine the heat kernel 
explicitly and ii) by applying the known relation between the residua of the product ({s) T{s) 



and B, 



n/2 



- — ^TTTTT = lim 
(47r)^/2 



s + 



n — d 



ns)as) 



n 



0,1,2,... . 



(4.13) 



(4.14) 



Both calculation schemes give the same values for the coefficients Bn/2 with 7^ : 

5i/2 = v^, fii = -2g, 53/2 = v^g^ B2 = -^^q\ 
B,/2 = ^q\ i?3 = -^g^ Br/2 = ^q^ .... 

Equations ()4.7j) - ()4.9|) give for the coefficient Bq 

Bo = 0, (4.15) 
while employment of the zeta function ()4.5p and relation ()4.13|) leads to the nonzero value 

1 



B. 



(4.16) 
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The different values of Bo in Eqs. fl4.15p and ()4.16p are the consequence of the noncompact 
configuration space, M^, in the problem at hand. In the case of compact manifolds the 
coefficient Bo is equal to the volume V of the manifold: Bo = V. Definition of Bq for a 
noncompact configuration space requires a special consideration jl^ |23| . The correct value 
of Bq in our case is formula (j4.15p derived from the explicit heat kernel (|4.7j) - (j4.9|) . The 
point is, the contribution of unbounded space is subtracted from the spectral density 
(14. 4p . Hence this contribution is subtracted also from the multiplier K^'^~^\t) in Eq. (|4.7p 
and, due to the factorization this formula, the contribution of is automatically subtracted 
from the total kernel K{t). 

When calculating the zeta function C^l-^) the spectral density ()4.4|) with subtracted 
M'^-contribution has been used also. However in Eq. ()4.5p . unlike Eq. ()4.7|) . there is no 
factorization of the contributions of waves propagating in directions parallel to the plasma 
plane and in normal to this plane direction. It is this point that lead to nonzero value of Bq 
in Eq. KT^ . 

B. TM- modes 

In the case of TM-modes there are two branches of spectrum: i) the photon branch (|3.34p 
and ii) the surface plasmon branch ()3.35|) . 

By making use of Eqs. ()3.27|) and ()4.4|) we obtain for the photon zeta function in the 
TM-sector of the model 



The integration in ()4.17|) can be easily done by introducing first the polar coordinates 




(4.17) 



k = rsm(p, p = rcosip, dk dp = r dr d(p, 
<ip < 7r/2, < r < oo. 



(4.18) 



In terms of new variables formula ()4.17|) acquires the form 




(4.19) 



Now it is natural to make the following substitution of the variables 



cos(p = x, —sm(pd(p 



dx . 



(4.20) 



As a result one gets 




(4.21) 



The integral over r in (j4.2H) exists in the region 



1/2 < Re s < 3/2 . 



(4.22) 
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and the subsequent integration over x can be done for 

1/2 < Re s < 1. 



(4.23) 



Thus the integral representation ()4.2H) can be analytically extended all over the complex s 
plane. This yields 



-,2-2s 



TM/ 
ph \ 



[s-m-s] 



r(3/2-s)rGs-i/2) 



^vr (1 — s)(2 — s) cosvrs 



. (4.24) 



Let us proceed to the construction of the heat kernel for the photon branch of the spectrum 
in the TM-sector of the model at hand. By making use of Eqs. dHSll), (jOl) , (031) , and (lOjl . 
we obtain 



27r2 



kdke 



dpe 



-p-'t 



(4.25) 



'0 JO 

Further we again introduce the polar coordinates ()4.18|) and then make the change of vari- 
ables ()4.20j) . Finally calculations are reduced to the substitution of the term r"^'' in Eq. 
(jOT|l by 



where 



CO ^2^-rH 



27r2 



dx (2x 



1)1 x) 



(4.26) 







+ q^x"^ 



dr 



In Tiqx ^2,j.2^ 

2 V T ~ 



1 — erf ( qx\/t 



1 /tF vrgx 
2V i ~ ~2~ 



2, 2 ^ 2\qx^tf'^+^ 

(2k 
^ fc=0 ^ 



1^ 



(4.27) 



We have used here the series representation for the function erf (qxy/t) j2J, UL 

Equations ()4.13p . ()4.24|) . on the one hand, and the explicit form of the heat kernel, Eqs. 
(I4.26|) and ()4.27|) . on the other hand give the same value for the coefficients i?„/2 in expansion 
()4.12|) for all n save for n = 0, 1: 



5i = -|g, 53/2 = 0, 52-^-3 
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q , 



5/2 



^9 , ^3 



8 „5 R , _ /,6 



(4.28) 



The explicit form of the integrated heat kernel in question ()4.26|) and ()4.27|) shows that two 
first coefficients in expansion ()4.12|) vanish while formulas ()4.26|) and ()4.27|) give 



-Bo = 5 Bi/2 

q 



TT , 



(4.29) 



This discrepancy implies that calculation of the first heat kernel coefficients through the 
respective spectral zeta function requires a special consideration in the case of unbounded 
configuration space. 

Let us address now the spectral functions for the spectrum branch ()3.35p generated by 
the surface plasmon. Equations ()3.35|) . ()4.1|) . and ()4.3p give 



TM/ 
sp \ 



d'k <^(k) 



(27r)2 



k dk 



(4.30) 
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The convergence of this integral in the region — requires Re s < 1, but when A; — oo it 
exists only if Re s > 2. Thus for the spectrum ()3.35p it is impossible to construct the zeta 
function by making use of the analytical continuation method, since one cannot define this 
function in any finite domain of the complex plane s. 

However it turns out that for this branch of the spectrum the heat kernel can be con- 
structed explicitly. Indeed, by making use of Eqs. (|4.2|) . (|4.3j) . and (j3.35p we obtain 



sp 



exp 



sp 



(k) 



k dk 



exp 



(4.31) 



The change of variables 



4A;2 



k dk = { ^ -\ — ] dx 
2 



reduces the integral ()4.31|) to the form 



K™ 



sp 



it) 



dx f X 1 
2^l? + 2'^ 



tx 



1 



+ 



1 



27iqH^ A-nt 



(4.32) 



The first term in Eq. ()4.32|l is absent in the standard expansion ()4.12|1 . and the second term 
in ()4.32|1 yields B1/2 = 2y^. The rest of coefficients -B„/2 with n 7^ 1 equal zero. Thus the 
surface plasmon with the spectrum ()3.35j) is a simple, and at the same time of a physical 
meaning, model that has no spectral zeta function, but the the respective (integrated) heat 
kernel exists though with a nonstandard asymptotic expansion. It is not surprising because 
we are dealing here with a singular point interaction described by the matching conditions 



V. LOCAL HEAT KERNEL 



In the TE-sector of the model under consideration one can construct, in an explicit form, 
the local heat kernel or Green function 



K{r,r';t) =< r|e-*^|r' >= ^ <(r)<^„(r')e^""* 



(5.1) 



We are using here the same notations as in Eq. (j4.2|) and V9„(r) are normalized eigenfunctions 
in the spectral problems at hand 



Lipn{r') = Xn^n{r') L = -A 



(5.2) 



From the definition (j5.1|) it follows, in particular, that the local heat kernel obeys the heat 
conduction equations with respect the variables (r,t) and (r',t) 



Ar- 



di 
d_ 
di 



K{T,r';t) = 



(5.3) 
(5.4) 
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and initial condition 

K{r, r'; t) 6{r - r ), when t 0+. (5.5) 
Obviously in the case of TE-modes K{r, r'; t) can be represented in a factorized form 



K{r, r'; t) = k'^^-^'^ (s, s'; t) ■ K{z, z'; t), r = (s, , 



where K^l^ ^■*(s, s'; t) is the free heat kernel in directions parallel to the plane z = 0: 



(d=2) , 



s, s'; t) = exp 

' ^ Ant ^ 



it 



(5.6) 



(5.7) 



For constructing the heat kernel K{z, z'; t) along the infinite 2-axes we shall use the tech- 
nique of integral equations developed in our previous paper [28^. These equations naturally 
arise when representing the Green function to be found in terms of heat potentials of simple 
and double layers. With respect of its first argument z the heat kernel K{z, z'; t) obeys the 
matching conditions ()3.13|) and when t — > 0+ it satisfies the initial condition ()5.5p . 

In what follows it is convenient to represent the heat kernel we are looking for in terms 
of four components depending on the range of its arguments 



K{z,z';t) = < 



(5.8) 



' K_+{z,z';t), z <0, z' >0, 

K++{z,z'-t), z,z'>0, 

K+_{z,z';t), z>0, z' <0, 

^K__{z,z';t), z,z'<0. 

We represent the heat kernel K{z, z'\ t) in terms of heat potentials of a simple layers [2 

K.+{z,z';t) = [ dTKo{z,0;t -T)ai{T;z') , z<0,z'>0, (5.9) 
Jo 



K+^{z,z';t) = Ko{z,z';t)+ / dr Ko{z,0;t — t) a2{T; z') , z, z' > . 

Jo 



(5.10) 



where ai^r; z') and a2{T; z') are the densities of the heat potentials to be found and 
Ko{z,z']t) is the free heat potential (Green's function) on an infinite line 



Koiz,z';t) 



exp 



[z — z 



/\2 



At 



(5.11) 



Substituting Eqs. ()5.9|) and ()5.10|) in matching conditions ()3.13p we obtain 

/ rfriro(0,0;t-r)[ai(r;z') -«2(r;/)] = ^o(0,z';t), ^' > , (5.12) 
Jo 

2^^(z = 0, z'; t)-Aqf drKoiO, 0; t - r)ai(r; z') = aAr; z') + a2(r; z'), z' > 0.(5.13) 
dz Jo 

When deriving Eq. (j5.13|) we have taken into account that the derivative of the single layer 
potential has a jump at the interface, namely, let V{z] t) is the heat potential of a single 
layer with the density z^(t): 



Viz;t) 



dTKo{z,0;t-T)u{T) 



(5.14) 
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then one can easily show that 2^ 2^ 



^{z = 0-;t) = --Ht), -(. = 0-;t) = -.(t). 



(5.15) 



On substituting in Eq. ()5.12j) the exphcit form of the free heat kernel ()5.1H) we arrive at 
the Abel integral equation the exact solution of which is known (3fll l3ll | or it can be derived 
by the Laplace transform 



z I z 

a\(t\ z) — a2{t; z) = — exp I ) , z > . 

W Ant \ 4t 



(5.16) 



In what follows it is convenient to apply to Eqs. (j5.1Hj) and (j5.1f)j) the Laplace transform 

m= / e-^'f{t)dt. (5.17) 
Jo 

As a result we obtain the set of two linear algebraic equations for the functions ai{p; z') and 
02 (p; z'): 



ai{p; z') - a2{p; z') = e ""'^^ / > 

2g ^ 



ai{p] z) + a2{p:, z) = e ""'^ - ^ai{p; z'), z' > . 



From here we deduce 



ai{p; z!) 



q + Vp 



«2(p; z') 



q + Vp 



-"'v^, z' >o. 



(5.18) 
(5.19) 

(5.20) 



By making use of the Laplace transform ()5.17|) the initial formulas ()5.9p and ()5.1U|) can be 
rewritten in the form 



K^+iz,z';p) = -— e-l^lv^«i(p;z'), z<0, z'>0, 

^Vp 

K^+{z,z';p) = ^e-l^-^'lv^+-^e-^v^a2(p;^'), z, z' >() . 

^Vp ^Vp 

Substitution of Eq. (PT:^ into Eqs. (IHTTTD and dlT^ yields 

1 



K-+{z,z';p) 
K++{z,z';p) 



2(g + v^) 
1 



<\z\+z')Vp z<0, z' > 



(5.21) 
(5.22) 

(5.23) 



2^^^ '2^{q + y/p) 

The inverse Laplace transform [i^ S^] applied to Eqs. ()5.23p and ()5.24|) gives 

\z\ + z' 



z, z' >0. (5.24) 



K_+{z,z'-t) = ifo(kl +^',0;t) - le-^^l^l+^'^+^'^erfc (^gv^ + 
;z < 0, / > , 



2y/i 



(5.25) 



K++{z,z';t) 



Koiz, z'; t) - I e«(-+^')+'?'* erfc (^qVi + , z, z' > . (5.26) 
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By making use of the same technique one can derive for the components Kj^-{z, z' \t) 

and K {z,z']t) the formulas which are analogous to Eqs. ()5.25|) and ()5.26p . It turns out 

that all four components of the heat kernel K{z^ z'\ t) (see Eq. ()5.8p ) can be represented in 
a unique way, namely: 

K{z, z'- 1) = Koiz, z'- 1) - I e'^i\A+\A)Wt ^^^^ ^ ^^^T^) ' 

-oo < < cx), zi^ 0, z' ^ 0. (5.27) 



In literature |33l . I34| the fundamental solutions to the heat conduction equation or to the 
Schrodinger equation with the Dirac 5 potential are represented, as a rule, in an integral 
form 

/■oo 

K{z, z] t) = Koiz, z] t)-q e""" Ko{\z\ + \z'\ + u,0;t) du . (5.28) 

Jo 

It is worth noting that the formulas (|5.27|) and (|5.28p have the sense also at the points z = 
and z' = 0. 

Thus the technique of integral equations developed on the basis of heat potentials in 
Ref. jl^ affords a regular and simple method for deriving, in an explicit form, the heat kernel 
in the problem under consideration. In Ref. |3^ the representation ()5.27p was obtained in 
the framework of another approach to this problem. 

Proceeding from Eq. (j5.28j) one can easily obtain the traced heat kernel (j4.9|) by evaluating 
the integral 

/oo 
dz [K{z,z;t) - Ko{z,z;t)] 

/OO POO 
dz due~''''Ko{2\z\+u,0;t). (5.29) 
■oo Jo 

Substitution to this equation of the exact form of the free heat kernel ()5.11|) gives 

(2z + u)^' 



Kit) = 1= / dz [ du e-«" exp 

V vrt Jo Jo 



At 



Thus, in order to find the integrated heat kernel K{t) one has to calculate the Laplace trans- 
form of the function erfc(M/ (2v^)). By making use of the table of the Laplace transform 
one finds 

i^(t) = le*«%rfc(gVt)-^. (5.31) 

The first term on the right-hand side of Eq. ()5.3H) exactly reproduces Eq. ()4.9|) while the 
constant term (-1/2) should be removed by introducing appropriate renormalization when 
calculating the trace over the continuous variables. 

In the TM-sector of the model under study the construction of the local heat kernel 
proves to be much more complicated. First of all, the factorization equation ()5.(j|l does not 
hold anymore. However, the technique of integral equations ^| can be applied here again, 
but the pertinent integral equations will be defined on the whole plane z = 0, i.e., they 
are two-dimensional integral equations. In order to go along this line, preliminary one has 
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to remove the dependence on the spectral parameter (a;(p, k)/c)^ from the first matching 
condition ()3.12|) . For the heat kernel K{r, r; t) this condition can be transformed as follows 



2g^(s,^ = 0,r';t), (5.32) 



where the notations r = (s, 2;), s = {x,y), and ()3.14|1 are used as before. 

The validity of this relation can be checked easily if one takes into account that the 
eigenfunctions (p{r) e~^"* entering the definition of the heat kernel (jS.lj) obey Eqs. (j5.3|) and 

-A (^„(r) e-^"* = A„ (/.„(r) e"^"* = -|- ^„(r) e"^"* . (5.33) 

ot 

Here the eigenvalues A„ are equal to (ti;(p, k)/c)^. The matching condition (j5.32|) involves 
the space and time derivatives, i.e., so-called the skew derivative. 

The Green function K{r, r'; t) obeys the heat conduction equations ()5.3|) and ()5.4|1 . hence 
the matching condition ()5.32|) can be obviously rewritten in the form 

OK 

[ArK{s, z = 0, r'; t)] = 2q — {s,z = 0, r'; t) . (5.34) 

An interesting spectral problem arises here when we confine ourself to the one-dimensional 
problem with {uj{p)/cY = p^'- 

— r Lfiz) = (f(z) , — oo<2;<oo, 2:7^0, (5.35) 

dz 

[^\z = 0)] = 0, (5.36) 

[ip"{z = 0)] = 2q^'{z = 0). (5.37) 

At first sight, this spectral problem is completely different from the analogous one for the 
6 potential (see Eqs. 1)3.111) . ITT^j) . and I3.13|) ). However the respective eigenfunctions (Pp{z) 
and ^p{z) are connected by relation 

^;(z) = $,(z) . (5.38) 

Both the problems have the same positive continuous spectrum 

0</<oo. (5.39) 

Furthermore the respective phase shifts, and consequently the scattering matrices, coincide 
(see Eq. ()3.21|) and Eq. ()3.23p with k = 0). From here we infer immediately that these 
spectral problems have the same integrated heat kernels defined by Eq. ()4.9j) . while the local 
heat kernels are obviously different. 

The local heat kernel for the spectral problem ()5.35|) - ()5.37|) can be derived in the same 
way as it has been done above in the spectral problem with 6 potential by making use of 
the integral equations. 

As far as we know, a couple of spectral problems possessing such interesting features are 
found for the first time. 
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VI. CONCLUSION 



The plasma sheet model investigated here proves to be very interesting and instructive, 
first of all, from the stand point of spectral analysis. The spectrum of the model contains 
both continuous branches and bound states (surface plasmon). It is remarkable that for 
the latter the spectral zeta function cannot be constructed at least by the standard analytic 
continuation method. At the same time the integrated heat kernel is found in an explicit form 
for all branches of spectrum. On the whole this heat kernel has the asymptotic expansion 
of a noncanonical form due to the singular point interaction. 

By making use of the technique of integral equations developed by us earlier ji^ the local 
heat kernel in the TE-sector of the model is found in an explicit form. By the way, here 
a new method is demonstrated for deriving the fundamental solution to the heat conduc- 
tion equation(or to the Schrodinger equation) on an infinite line with the 5-like source. In 
principle the integral equations technique is also applicable to construction of the local heat 
kernel in the TM-sector of the model. 

For the heat equation on an infinite line with the (5-source a nontrivial counterpart is 
found, namely, a spectral problem with point interaction, that possesses the same integrated 
heat kernel. However the local heat kernels in these spectral problems are different. The new 
boundary problem is not selfadjoint, the matching conditions at the singular point can be 
represent ed in two forms: by equation containing the spectral parameter or ii) by equation 
with second spatial derivative (or the first time derivative). As far as we know such a couple 
of spectral problems is found at the first time. 

Thus the spectral analysis of the model under consideration is accomplished in full. 

Without question, the spectral analysis of the same model describing the plasma layer of 
other forms, for example, circular infinite cylinder or sphere, is also of interest. 
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